Magnetic properties of a metal-organic antiferromagnet on a distorted honeycomb 

lattice 
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For temperatures T well above the ordering temperature T„ = 3.0±0.2K the magnetic properties 
of the metal-organic material Mn[Ci H 6 (OH)(COO)]2x2H 2 O built from Mn 2+ ions and 3- hydroxy- 2- 
naphthoic anions can be described by a S = 5/2 quantum antiferromagnet on a distorted honeycomb 
lattice with two different nearest neighbor exchange couplings J 2 ~ 2Ji w 1.8 K. Measurements of 
the magnetization M(H, T) as a function of a uniform external field H and of the uniform zero field 
susceptibility x(T) are explained within the framework of a modified spin-wave approach which 
takes into account the absence of a spontaneous staggered magnetization at finite temperatures. 
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I. INTRODUCTION 

In recent years the role of fluctuations, spatial 
anisotropy and frustration in low dimensional quantum 
magnets has been intensely studied, both experimentally 
and theoretically^ For a comparison of experiments with 
theory it is crucial to have well defined crystalline mate- 
rials where one or several parameters can be varied ex- 
ternally in order to obtain quantitative predictions for 
physical observables. Moreover, in order to observe in- 
teresting magnetic many-body effects it is essential to 
have materials where the magnetic moments are cou- 
pled via sufficiently strong exchange interactions. These 
conditions are met by transition metal oxides such as 
cuprates, vanadates, copper-germanates, or manganites, 
which have been the subject of many works. However, in 
these materials it is rather difficult to control externally 
microscopic parameters such as the precise values of the 
exchange interactions or the lattice geometry by chang- 
ing the chemical composition in a well defined manner. 
This problem tends to be less severe in magnets based on 
metal-organic materials, which offer more possibilities of 
modifying some constituents chemically and thereby tun- 
ing the properties by a crystal engineering strategy. The 
challenge is then to find metal-organic magnets where 
the magnetic moments are coupled sufficiently strongly 
to exhibit interesting collective effects. 

These effects arc of particular importance in low- 
dimensional magnets, e.g. 2D layer structures with strong 
magnetic couplings within the layers and weak interac- 




FIG. 1: View along the 6-axis of the metal-organic quantum 
magnet Mn[Ci H 6 (OH)(COO)] 2 x2H a O. Bold lines show ex- 
change paths Mn — O — C — O — Mn. The unit cell, denoted by 
the parallelogram, contains four crystallographically equiva- 



tions between the layers. Such layer structures can be 
built up chemically from spin-bearing metal ions, which 
are connected by short bridges, being separated by or- 
ganic fragments of considerable size, see Fig. ^ Moti- 
vated by these considerations we synthesized transition 
metal complexes of o-hydroxy-naphthoic acids. The crys- 
tal structure of Mn[CioH 6 (OH)(COO)] 2 x2H 2 (system- 



FIG. 2: Chemical formula of Mn[Ci H6(OH)(COO)]2x2H 2 O. 



atic name: manganese(II) 3-hydroxy-2-naphthoate dihy- 
drate, Fig. is of particular interest, because the Mn 2+ 
ions form a distorted honeycomb lattice (Fig. |3J). For 
the redetermination of the crystal structure, pale brown 
crystals were slowly grown by diffusion of an aqueous so- 
lution of Na[Ci H 6 (OH)(COO)] into an aqueous MnS0 4 
solution with a buffer layer of water. The single crys- 
tal X-ray analysis confirmed the previously determined 
structure^ with a higher precision. The compound crys- 
tallizes in the monoclinic space group P2\/c with the 
lattice parameters a = 17.191(4) A, b = 7.3448(10)1, 

c = 15.5279(17)1, (3 = 101.964(8)°, V = 1918.1(5) I 3 . 3 
The unit cell contains four crystallographically equivalent 
Mn 2+ ions. 

The coupling layer, parallel to the (be) plane, contains 
the Mn 2+ ions, the COO~ and OH groups as well as wa- 
ter molecules. The isolating layer, having a thickness of 
about 12 1 consists of the organic naphthalene moieties. 
These naphthalene moieties are only bound together by 
van der Waals contacts between C and H atoms. The 
relative weakness of these interactions is reflected by the 
morphology of the crystals: the crystals grow in (6) and 
(c) direction much faster than in (a) direction, thus form- 
ing thin plates parallel to the (be) plane. 

The magnetism is due to the S = 5/2 manganese ions 
which form a distorted honeycomb pattern parallel to the 
(6c) planes. Neighboring ions are connected by carboxylic 
groups, which provide an Mn — O — C — O — Mn magnetic 
exchange path. There are two different exchange paths: 
the first path contains a single O — C — O unit, displayed 
in green in Fig. EI I n the second path (marked with blue 
color) the Mn 2+ ions are connected by two O — C — O 
moieties simultaneously. The honeycomb layers are well 
separated from each other; the closest distances between 
Mn 2+ ions of different layers are as large as 16.2821. 

The structure in Fig. [31 suggests that the magnetic 
properties of the material can be modeled by a spin 
S = 5/2 Heisenberg magnet on the distorted honey- 
comb lattice shown in Fig. 0] The exchange integrals 
J u = J(ri, Ti + S v ), v = 1,2, 3, couple the spin at a given 
site Ti to its nearest neighbors at rj + <5„. All exchange 
integrals J„ turn out to be positive, and |<5i| = = 
<5i = 5.131(4)1 and J\ — J 3 , due to the crystal sym- 
metry. A closer look at the crystal structure in Fig. [3| 



FIG. 3: View on the (6c) plane of the metal-organic quantum 
magnet Mn[CioH 6 (OH)(COO)] 2 x2H 2 Q. 



and a comparison with the distorted honeycomb lattice 
in Fig. 0| reveals that J2 acts along two exchange paths 
while J 1 results from a single exchange path. Therefore 
we expect J2 to be roughly twice as large as J\ . Further- 
more, the honeycomb lattice is bipartite, i.e., it can be 
divided into two sublattices, labeled A and B, such that 
the nearest neighbors of all sites belonging to sublattice A 
are located on sublattice B. Thus, for positive J v the sys- 
tem is not frustrated, and when quantum fluctuations are 
neglected the ground state shows classical antiferromag- 
netic Neel order. More generally, we expect long-range 
antifcrromagnetic order to persist in the quantum me- 
chanical ground state. Therefore, it should be possible 
to calculate the magnetic properties of the system within 
the usual spin-wave expansion, at least for temperature 
T = 0. Note that the actual structure shown in Fig. |3 
has an additional distortion in the x-direction, resulting 
in a primitive cell with doubled volume. Due to the low 
symmetry of the lattice the Dzyaloshinskii-Moriya inter- 
action might play an important role. However, we expect 
the corresponding energy scale to be small in compari- 
son with Ji and J2, so that in the first approximation 
we can neglect this effect. In the following we therefore 
always work with the magnetically equivalent Bravais lat- 
tice shown in Fig. 0] 

Measurements of the magnetization M(H, T) in a mag- 
netic field H are performed at finite temperatures T, 
where long-range antiferromagnetic order is ruled out 
by the Hohenberg-Mermin- Wagner theorem^ In this 
case the theoretical justification for the spin-wave ex- 
pansion in two dimensions is more subtle. As long as 
there is long-range antiferromagnetic order at T = 0, it 
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O sublattice A 
sublattice B 



FIG. 4: Distorted honeycomb lattice. The interactions be- 
tween spins are displayed as solid lines. The underlying mag- 
netic sublattice is a Bravais lattice and its primitive cell can 
be chosen to be the dashed parallelogram. The correspond- 
ing primitive vectors are a\ = ai& x and a 2 = a^cos^e^ + 
a 2 sin if e y . 



is reasonable to expect that the low-energy and long- 
wavelength physics is still dominated by rcnormalized 
spin- waves. The magnetic properties of square lattice an- 
tiferromagnets in the absence of uniform external fields 
have been thoroughly investigated in a classical work by 
Chakravarty, Halperin, and Nelson^ Less is known about 
the low-energy physics of two-dimensional quantum an- 
tiferromagnets subject to a uniform external magnetic 
field. The external field breaks the rotational symme- 
try of the Heisenberg antiferromagnet to 0(2), similar to 
the effect of an XY anisotropy in the XXZ model e 6 How- 
ever, the classical ground states of the two models differ 
substantially: whereas the XXZ model has a collinear 
ground state, an uniform magnetic field in a Heisenberg 
antiferromagnet leads to a canted classical spin configu- 
ration shown in Fig.JSJ The zero-temperature magnetiza- 
tion curve M(H, 0) of the square lattice antiferromagnet 
has been calculated a few years ago by Zhitomirsky and 
Nikun: 7 within the spin-wave expansion. For finite tem- 
peratures, M{H,T) has been extrapolated from numeri- 
cal diagonalizations of finite clusters. 8 We are not aware 
of any analytical calculations in the literature of M(H, T) 
for two-dimensional quantum Heisenberg antiferromag- 
nets at T > 0. In this work, we calculate M(H, T) using a 
modified spin-wave approach™ which takes the absence of 
a spontaneous staggered magnetization at finite tempera- 
tures into account. Our theoretical results for the magne- 
tization curves as well as for the zero-field susceptibility 
X(T) show a satisfactory agreement with our measure- 
ments for the compound Mn[Ci H 6 (OH)(COO)] 2 x2H 2 O. 

The rest of the paper is organized as follows. In Scc.lTTl 
we review the formalism of the spin- wave expansion for 
non-collinear spin configurations. In Scc. lIIII this method 
is applied to an antiferromagnet on a bipartite lattice in 
the presence of a uniform magnetic field. Expressions 
for the magnetization, the staggered magnetization and 




FIG. 5: Spin configuration (Si) in the classical ground state 
of a two-sublattice antiferromagnet . The dashed arrows rep- 
resent a uniform magnetic field Be x in the ^-direction and a 
staggered magnetic field (iB s e z in the z-direction. The small 
solid arrows represent the vectors of a "co-moving" basis that 
matches the direction defined by the local magnetization (Si). 
Not shown are the basis vectors e A = ejj = e y which point 
into the plane of the paper. 

the uniform susceptibility for the material of interest are 
obtained. We explain how a self-consistently determined 
staggered field is used to regularize divergencies at finite 
temperature. In Sec. II VI we present results and compare 
with experimental measurements. Finally, in Sec. we 
present our conclusion. 



II. SPIN WAVES IN NON-COLLINEAR SPIN 
CONFIGURATIONS 

In the presence of a homogeneous magnetic field an an- 
tiferromagnet on a bipartite lattice has a non-collinear, 
canted spin configuration as shown in Fig. [5j We choose 
a coordinate system such that the uniform external field 
points along the x-axis, and the staggered magnetiza- 
tion is directed along the z-axis. The low temperature 
physics is dominated by spin-wave excitations. To obtain 
their spectrum we should quantize the spin-operators in 
a spatially-dependent ( "co-moving" ) coordinate system 
that matches for each site the axis defined by the expec- 
tation value (Si) of the spin operator. 

More generally, the problem of calculating the spin ex- 
citations of a Heisenberg magnet subject to an arbitrary 
inhomogeneous magnetic field Bi can be formulated and 
solved in a coordinate-free vector notation^ Consider 
the general Heisenberg hamiltonian 

H = ^ JjjSj ■ Sj - g(j, B J~] Bi ■ Si , (2.1) 

i,j i 

where = J(rj,r 3 -) are some arbitrary exchange cou- 
plings, the sums are over all sites of a D-dimensional 
lattice consisting of N sites, and the Si are spin-S* op- 
erators normalized such that Sf — S(S +1). The last 
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term represents the Zeeman energy, where g is the gy- 
romagnetic factor and «b is the Bohr magneton. We 
assume that the external magnetic field is sufficiently 
strong to induce permanent magnetic dipolc moments 
rrii = gfj,-Q(Si), where (. . .) denotes the thermal equilib- 
rium average. It is then convenient to decompose the 
spin operators as Si — s}rhi + S-~, where ■ rhi = 0, 
and rhi = mi/\trii\ is a unit vector in the direction of 
mi. Substituting this decomposition into Eq. 1|2.1|) wc 
obtain H = + H x + H', with 



m = 

H 1 - = 
H' = 



\Y,JijSi-Sf, (2.3) 



^^•(^-^J^I'm,), 



(2.4) 



where hi = g^Bi. Note that H' describes the coupling 
between the transverse and the longitudinal spin fluctu- 
ations. The classical ground state energy E$ is obtained 
by replacing S- — > S in Eq. (|2.2|l and by finding the 
configuration {rhi} that minimizes the resulting classi- 
cal hamiltonian 



H c 



S 2 



Sj^hi 



(2-5) 



A necessary condition for an extremum of Eq. (|2.5(1 , tak- 
ing into account the constraints rh 2 = 1, isiS 



h; 



Sy ./, ; mj (i. 



3 



(2.6) 



For given hi and , this is a system of non-linear equa- 
tions for the spin directions rhi in the classical ground 
state. Using Eq. (|2.(i|) . the part H' of the hamiltonian 
describing the coupling between transverse and longitu- 
dinal fluctuations can be written as 



H' 



,)(s-s). 



(2-7) 



1,3 



Let us expand the transverse components of S± in a 
spherical basis, = \ Ylp=± S^ v e?i , with the spherical 
basis vectors eP i = e\ + ipe 2 , p = ±, where {e 1 ,^ 2 ,^} 
is a local orthogonal triad of unit vectors. The transverse 
part of our spin hamiltonian can then be written as 



(2.8) 



The basis vectors e},e 2 are not unique: any rotation 
around rhi yields an equally acceptable transverse basis. 

So far, no approximation has been made. To obtain 
the spin- wave spectrum, we expand the spin operators in 



terms of canonical boson operators bi as usua l 11 ' 12 , Sf = 
S-b\b l and S+ = (Sr)t = ^256. (l + O^- 1 )). Within 
the linear spin-wave approximation the hamiltonian be- 
comes H ps Eq 1 + H2, where Eq 1 is the minimum of the 
classical hamiltonian H cl in Eq. (|2.5() . and 

h, = fE J a t( e « + ■ e 7) h h + k ■ 4) bfa 

+ (et ■ e+) b\b] + (er ■ ej) b,b t ] 

g 



{hi ■ rhi)b\bi 



(2.9) 



Note that the contribution from H' in Eq. I|2.7J) is of order 
S 1 / 2 and hence can be neglected within linear spin-wave 
theory. Eq. I|2.9|l together with Eqs. 1|2.5() and (|2.6|l com- 
pletely determine the spin- wave spectrum of any Heisen- 
berg magnet in an arbitrary inhomogeneous field. 



III. SPIN WAVES IN THE DISTORTED 
HONEYCOMB LATTICE 

A. Classical ground state 

Let us apply the general formalism of the previous sec- 
tion to our bipartite lattice antifcrromagnct in a uniform 
external magnetic field Be x along the x axis. We denote 
by e a fixed unit vectors in direction a = x,y, z. For 
technical reasons we introduce an additional staggered 
magnetic field QB s e z in the z-direction, where Ci = 1 
if tv; belongs to the A-sublattice and & = — 1 if r; be- 
longs to the B-sublattice. This auxiliary staggered field 
will be determined self-consistently in Sec. lIII Ol to insure 
a vanishing staggered magnetization at finite tempera- 
tures. The total magnetic field is thus 



hi = gfiB[Be x + &B 8 e z ] 



(3.1) 



The classical ground state configuration is then rhi = 
(i cos 8e z + sin 8e x , as shown in Fig. [5] 

For convenience we introduce the notation no = cos 9 
and mo = sin 6*. Physically, mo corresponds to the clas- 
sical limit (S — > 00) of the normalized uniform magneti- 
zation 



1 

~NS 



S t ) 



(3.2) 



while no corresponds to the S 
ized staggered magnetization 



00 limit of the normal- 



1 

~NS 



Eg 



Si) 



(3.3) 



By symmetry, the uniform magnetization points into the 
x-direction, while the staggered magnetization points 
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into the z-direction. The natural dimensionless mea- 
sure for the strength of the fields is h — xoffMB-B and 
h a = Xo9^bB s , where \o = (2Jo5 , ) _1 is the classical 
uniform susceptibility. Here Jo = ^2„Jv is the k = 
component of the Fourier transform of the exchange cou- 
plings 

j k =j2 e ~ ikSuJ v ( 3 - 4 ) 

For the special choice of the field hi given in Eq. (|3.1|) 
our general Eq. I|2.6|l reduces to the simple relation 

h = mo[l + h s /no] , (3.5) 

which together with n\ + — 1 determines the classical 
Neel order parameter no and the classical uniform mag- 
netization nio as functions of the fields h and h s . Note 
that rh.A • ttib = — n 0i an( l w ith the special transverse 
basis shown in Fig. 

e A • e B = 2[5 PiP ml + dp-p/ml 1 ] , (3.6) 
rh-A • = 2ipnomo = — rhe • e A . (3.7) 



ky 



b 2 



FIG. 6: Reduced Brillouin zone of the distorted honeycomb 
lattice. The primitive vectors are b\ = ^ (sin tp e x — 

cos <pe y ) and 62 = - where ai, a2 and the angle 95 
are defined in Fig. 2 



B. Spin- wave dispersion 

To obtain the spin-wave dispersion, we must diagonal- 
ize H2 in Eq. I|2.9(l for the special ground-state spin con- 
figuration discussed above. Therefore, we first perform 
Fourier transformations separately on each sublattice, 

k = v/-|xy fc ' riafe - f ° rr ' eA ' ( 3 - 8a ) 

k 



\ A Ah- 



N ^ 

k 



r >b k , for n £ B , (3.8b) 



where the wave- vector sums are over the reduced (mag- 
netic) Brillouin zone of the honeycomb lattice shown in 
Fig. El With the above definitions we obtain 



H 2 = JoSY,[Ma{a k + b{b k )+B k b_ k a k +B* k alb i _ k 



+C k b{a k + C* k alb k ] 



(3.9) 



where A = 1 + 2h s /no, B k = ii^J k /Jo, and C k = 
iTiQj k /Jo. On a honeycomb lattice J k — | Jk|e** fc is com- 
plex, so that B k = \B k \e^» and C k = |C k |e^ fc . Using 
<f>- k = —(j) k , it is easy to see that these phase factors 
can be removed from Eq. I|3.9|l via the gauge transforma- 
tion a k = e l ^ h a k . Introducing then new canonical boson 
operators 



Cfccr = -k= [a k + ab k ] , a = ±1 



(3.10) 



the hamiltonian (|3.9|l assumes the block-diagonal form, 



H2 = ^J2[( A +^ c ^( c l c ^+'-^-^' 



ka 



The diagonalization is completed by means of the Bogo- 
liubov transformation, 



c k<J 

-fco 



u ka -ov ko 
-av ka u ka 



— ko 



where 



u ko 



Vka = 



lA + (j\C k \ +e ka 



2e ka 



! A + a\C k \ -e ka 



2e ko 



with the dimensionless energy dispersion 

e ka = y/(A + <r\C k \) 2 - \B k \< 
Defining j k = J k /Jo, we may write 



(3.12) 

(3.13a) 
(3.13b) 

(3.14) 

1/2 



'fc~= (1 + — +cr|7fe|)(l + — -v{nl-rnl)\"i k \) 
L no n J 

(3-15) 

In terms of the new operators d ka the quadratic spin 
wave hamiltonian H2 is diagonal, 

H 2 = JoSY,Ukad ] ka d ka + \ [tk* -(A + a\C k \)} \ . 

(3-16) 

The low temperature properties of the magnet are deter- 
mined by the long-wavelength behavior of the spin-wave 
dispersions, which follow from the expansion for small fc, 



I7fc| 



1 _ \ 12 k » A «P k P > (3-17) 



mi 
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where A is a matrix with elements 

A a p, = E ir(8u-eg)(8u'ep)—y y j (S v -e a ){S V '-ep) . 
Jo 



v.v J o 



(3.18) 

Since A is symmetric, an orthogonal basis can always be 
chosen such that A is diagonal, with eigenvalues A a . In 
this basis 



17*1 



1 E A a k 2 a . 



The matrix A is positive, since 



l7fc|<E 



(3.19) 



(3.20) 



velocities c Q 
are thus 



JoS£ a along the two principal directions 



c / J\ Jo 
= S^j — a u 

= S\JlJ\J2 a,2 sin ip . 



(3.21) 
(3.22) 



Note that for J 2 — > the velocity c y vanishes, so that 
the system becomes one-dimensional, as is obvious from 
Fig. On the other hand, for J x — > both velocities 
vanish, because in this limit the system consists of de- 
coupled dimers. 

For h s — only the mode is gapless for k — > 0, 
while the mode Ef,+ has the gap 2mo- To give a more ex- 
plicit form for the long- wavelength spin- wave dispersions, 
we further assume h s -C uq. Then 



n 



n 



S ' (Iqfcq) 



1/2 



4m 



9 4/i s , 9n 
n 



-{n 2 - 2m 2 )^(i a k a ) 



1/2 



(3.23) 



(3.24) 



For no — > the expansion H3.23JI is not appropriate any 
longer and for /i s = the dispersion becomes purely 
quadratic at no = 0. Before this happens, there is a criti- 
cal field < h* < 1 at which the curvature of the disper- 
sion efc_ changes sign. The positive curvature for h > h* 
results in an instability of magnons towards a sponta- 
neous decay into two magnon states^ Furthermore, if 



an anisotropic exchange is considered, the anisotropy gap 
A is strongly renormalized by magnon interactions^^ 
As the influence of these instabilities on the thermody- 
namic properties is unclear at the moment, they will not 
be further considered in this work. 



C. Uniform and staggered magnetization 

We now calculate the leading spin-wave corrections to 
the normalized uniform- and the staggered magnetization 
as defined in Eqs. I|3.2|) and (|3.3J) . A standard expansion 
in powers of 1/S gives 



where the last equality assumes that all couplings have 
the same sign. We can thus define effective length scales 
l a by setting A a — l 2 a . For a £>-dimensional hypercubic 
lattice with lattice spacing a we have (? a = a 2 /D. For 
our honeycomb lattice shown in Fig. 0] with |<$i| = | <5a | 
and Ji = J3 the eigenvectors of A are parallel to the x- 
axis and the y-axis, with corresponding eigenvalues I 2 , = 
{J\/2Jo)a\ and £ 2 = (2Ji J 2 / J§)a 2 sin 2 tp. The spin- wave an( j 



where 
F(h,h s ) 



Here n ko 



N 



T 
1 

"0 



E 



1 

N 



nh s F(h, h s ) 



1 — mpm 



S 

I(h, h s 



S 



(3.25) 
(3.26) 



47fe| 1 



E 

ho 



2h s 
no 



2K 
n 



(3.27) 



crm |7 fe 



= \ p JoSe h „/T 



(3.28) 

is the Bose function. 



The parameters n and m on the right-hand sides of 
Eqs. H3.25Kl3~2*5|) are determined as functions of the fields 
h and h s by Eq. I|3.5(l and rig + m§ = 1. Note that for 
S — > 00 the solutions of Eqs. 13.25J1 and H3.26(l correctly 
approach n — no and m = mo: in this limit Eq. 
reduces to Eq. (|3.5|l . while Eq. I|3.26|l simply becomes an- 
other way of writing n\ +nip = 1. In the thermodynamic 
limit, we transform Brillouin zone sums to integrals ac- 
cording to 



/ \ - N- 

k 



BZ 



d 2 k 



(3.29) 



where V u = aia2siny is the area of the magnetic unit 
cell in real space and the integral is over the reduced 
Brillouin zone shown in Fig. |SJ 

At T = and h s = expressions similar to (|3.25() 
and l|3.2(j|l have been discussed previously. 7 Only m(h) 
was given explicitly and a renormalization of the canting 
angle was found by considering spin- wave interactions. 
Yet, to a given order in 1/S it is easier to calculate m and 
n directly as derivatives of the free energy with respect to 
h and h s . Very recently, the renormalized canting angle 
was also used to analyze the behavior of n(h) at T = 
for a more complicated geometry^ 

At any finite temperature the integral I(h, 0) is in- 
frared divergent in two dimensions, signaling the absence 
of long-range antiferromagnetic order, in accordance with 
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the Hohenbcrg-Mermin- Wagner theorem. 4 At first sight, 
it thus seems that the finite-temperature magnetization 
curve cannot be calculated within our spin-wave ap- 
proach. Fortunately, there is a straightforward way to 
obtain an approximate expression for the magnetization 
even at finite T. The crucial observation is that if wc 
set n = in Eqs. I|3.25[l and l|3.26[) . these equations can 
be interpreted as a condition for the staggered field h s 
that is necessary to enforce a vanishing staggered mag- 
netization. The solution h s = h s (h) as a function of the 
uniform field h is not a physical external staggered field, 
but an internal effective field that is generated by strong 
fluctuations. In fact, the field h s (h) is nothing but the 
Lagrange multiplier introduced in Takahashi's modified 
spin-wave theoryAll It is well known that the internal 
field is related to a finite correlation length £, as we will 
further discuss in Sec. IIV CI Numerically, we calculate 
the uniform magnetization m(h, T) at finite temperature 
T by adjusting h s for fixed external field h such that the 
condition n = is ful filled in Eqs. and (|3~2BTl . Us- 

ing this h s (h) in Eq. (|3.25|1 then directly yields m(h,T). 

We must keep in mind that the staggered field h s does 
not respect the rotational symmetry of the original hamil- 
tonian, which for h = corresponds to a global 0(3) 
symmetry and for ft, > is reduced to a global 0(2) 
symmetry around the axis of the uniform field. With 
the parametrization that explicitly breaks this symmetry, 
we should therefore only calculate rotationally invariant 
quantities^ Below, we will find a disagreement between 
a rotationally invariant evaluation of the zero-field uni- 
form susceptibility and the slope ofdm/dh for h — > 0. We 
attribute this discrepancy to the fact that dm/dh\h^o 
does not respect the 0(3) symmetry in this limit. Gen- 
erally, we expect our approach for the finite temperature 
magnetization to be reasonable only for h > h s (h, T). In 
Sec. IIV Cl we will see that h s is exponentially small at low 
temperatures, such that h > h s (h,T) is fulfilled even for 
very small external fields. The condition h > h s (h,T) 
then roughly gives a limit of validity of our approach 
in terms of the temperature as T < 0.5JoS. The fact 
that the limits T — > and h — > do not commute 
in a modified spin-wave expansion was first noticed by 
Takahashi.il 



D. Uniform susceptibility 

In order to calculate the rotationally invariant uniform 
zero-field susceptibility per spin 

* = Yn ^ {Sl ' Sj) ' (3 - 30) 

i,j 

we set the uniform magnetic field B = in Ea. H3.1(l . In 
this case m = mo = and no = 1, so that we obtain a 
doubly degenerate mode in Eq. I|3.15|) with dispersion 

e k « = t k = v / (l + 2/i s ) 2 -|7fc| 2 , (3.31) 



and the expression for the staggered magnetization (|3.26() 
reduces to 




As explained in the previous section we use a self- 
consistently determined staggered field h s to enforce a 
vanishing order parameter n = 0. 

The susceptibility l|3.30|l can be written as 

X = j; (S q ,+ ■ S- q:+ ) q=0 , (3.33) 

where we have defined the linear combinations (a = ±) 

S qtt , = ±(S£ + <TS*) (3.34) 

of the Fourier-transformed spin operators on each sub- 
lattice 

S « /B = \Pti £ e~ iqri Si. (3.35) 

r,eA/B 

Next we decompose the susceptibility into a transverse 
and a longitudinal part 

X = X + -+X Z \ (3.36) 

where 

X + = 2T (^+^-9.+ + ^9,+^-9,+)q=o ' (3-37) 

X zz = f(S z q , + S z - q , + ) q=0 . (3.38) 

We map the spin operators H3.35(l onto canonical boson 
operators via a Dyson-Maleev transformationi 1 * 1 ^ and 
evaluate the thermal expectation values of the noninter- 
acting state using the Wick theorem. Then the trans- 
verse susceptibility H3.37(l is proportional to the right 
hand side of Eq. (|3.32(l , and thus vanishes if we require 
ri = 0. Therefore, in our approximation only the lon- 
gitudinal part contributes to the rotationally invariant 
uniform susceptibility, 

k 

Apart from a different normalization, this result has been 
obtained previously in Takahashi's approach,— We eval- 
uate Eq. H3.39|l numerically in the thermodynamic limit. 

IV. RESULTS 

A. Zero temperature uniform and staggered 
magnetization 

In Figs. [7] and |H] we show results for the uniform and 
staggered magnetization at zero temperature. In two spa- 
tial dimensions, there are no divergent contributions to 
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FIG. 7: Normalized uniform magnetization m(h) for T = 
and h s — 0. The solid line is the zero-temperature magne- 
tization curve for the honeycomb lattice with 5 = 1/2 and 
Ji = J2. For comparison we also show the corresponding 
curve for a square lattice and exact results for a linear anti- 
ferromagnetic chaini- However, the 5 = 1/2 chain is critical, 
so it is not surprising that it is poorly described by means 
of the spin-wave theory. Note that for h s = the classical 
equation 113. 511 is simply mo = h. 
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FIG. 8: Normalized staggered magnetization n(h) at T = 
for honeycomb (solid line) and square lattice with S = 1/2 
(dotted line) . The classical equation no = \/l — h 2 is plotted 
for comparison. We also show the curves for the anisotropic 
cases J1/J2 = 10, 100. 



for the square lattice {V u / M x l y =l) within linear spin- 
wave theory^ For our distorted honeycomb lattice, we 
have 



the integrals in Eqs. (|3.25|l and l|3.26[l . indicating true 
long range order. We can thus set h s = and conse- 
quently mo = h. As the deviations from the classical 
curves are rather small for S = 5/2, we present the curves 
for the extreme quantum case S = 1/2. 

The uniform magnetization shows a positive curvature 
for all < h < 1 and lies generally below the classical 
straight linei This tendency is stronger for the honey- 
comb lattice and is even more pronounced for anisotropic 
exchange couplings with J\ J2. The number of near- 
est neighbors 2 = 3 for the honeycomb lattice is lower 
than for the square lattice (2 = 4), and in the limit 
J2 <C J\ the system is almost one-dimensional. The 
observed tendency thus simply corresponds to increased 
quantum fluctuations in low dimensions. Beyond the sat- 
uration field h = 1 the ground state has full collinear 
ferromagnetic order. This state as well as single magnon 
excitations above it are easily shown to be exact eigen- 
states. As the single magnon states become gapless at 
exactly the classical value h = 1, the saturation field is 
not changed by quantum fluctuations or magnon inter- 
actions. The limit h — ► 1 is reached with infinite slope in 
m(h). The leading behavior is given by 



TO. 



1 



Vu Sh 
UJ y nS 



In (A6h) 



(4.1) 



where Sh — 1 — h. This logarithmic asymptotics was 
first discussed in the language of Bose condensation of 
magnons below the saturation field-fi and was later found 



Vu 

*-X* II 



\2J 1 + J 2 f 



(4.2) 



which diverges for J\ — > or J2 — * and thus exempli- 
fies the increasing deviations from the classical curve for 
strongly anisotropic exchange couplings. 

The staggered magnetization in Fig. [S] shows a non- 
monotonic dependence on the applied uniform field. For 
vanishing h the staggered magnetization decreases as we 
lower the effective dimensionality. An external field ap- 
parently suppresses quantum fluctuations and n(h) first 
increases with h before it reaches a maximum and then 
vanishes for h — > 1 with infinite slope. The asymptotic 
behavior is given by 



V u VSh 



In (45h) 



(4.3) 



Interestingly, the quantum corrections to the staggered 
magnetization are positive close to the saturation field 
and the spin-wave result therefore intersects the classical 
curve. In a quasi one-dimensional situation (J 2 <C Ji), 
quantum fluctuations are strong and the leading order 
spin-wave theory, when pushed to the limit of validity, 
predicts a quantum disordered phase for small uniform 
fields. 
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FIG. 9: Uniform magnetization m(h) for the honeycomb lat- 
tice with S = 5/2 and Ji = J2 for two values of T. 



B. Finite temperature magnetization and 
susceptibility 

Magnetic measurements were carried out on a sin- 
gle crystalline sample of Mn[Ci H 6 (OH)(COO)] 2 x2H 2 O 
with a mass of meoNA = 0.65 mg using a Quantum De- 
sign SQUID magnetometer MPMS-XL. Isothermal mag- 
netization runs at temperatures between 2 and 200 K and 
fields up to 5 T were performed as well as measurements 
of the susceptibility in the temperature range 2 — 300 K 
for a magnetic field of 0.05 — 2 T£ 

In Fig. [5] we show theoretical magnetization curves 
m(h) for the honeycomb lattice with S = 5/2 and 
Ji = J 2 at different temperatures T. For T <C JqS the 
magnetization is almost linear throughout the entire field 
range. At intermediate temperatures m(h) has an S-like 
shape with a positive curvature at small fields h that 
changes to a negative curvature with increasing h. Simi- 
lar low-temperature behavior of the magnetization curve 
has been observed in a quantum Monte Carlo study of the 
two-dimensional Heisenberg antiferromagnet on a square 
latticed 

It turns out that the magnetization as well as the sus- 
ceptibility are not very sensitive to the ratio J2/ J\ as long 
as Ji and J2 have the same order of magnitude. Thus, 
we cannot determine the precise value of J2/J1, but our 
fits are compatible with the assumption J2 « 13\. 

In Fig. 1101 we show experimental data and theoreti- 
cal fits for the normalized uniform magnetization m = 
M/(NS) at different temperatures. The magnetic field 
H = 2JgSh is given in Tesla. Surprisingly, all experimen- 
tal curves are almost straight lines, whereas from Fig. [5] 
we would expect an upward bend of m(h) at higher tem- 
peratures. Fits for T = 2K and different ratios J1/J2 
invariably give Jo ~ 4K. Hence we assume J2 = 2Ji 
and fit the theoretical curve to the experimental data at 




H[T] 

FIG. 10: Magnetization m(H) of Mn[Ci H 6 (OH)(COO)] 2 X 
2H2O up to field H = 5T. Experimental data are indicated 
by squares (2K), circles (8K) and triangles (20 K). Theoreti- 
cal magnetization curves for honeycomb lattice with S = 5/2 
and Ji — 2Ji = 1.95 K are denoted by lines. 




100 200 300 



T[K] 

FIG. 11: Susceptibility X (T) of Mn[Ci H 6 (OH)(COO)]2 X 
2H2O. Circles are experimental data in a field of 2T. The- 
oretical fit for honeycomb lattice with J2 = 2Ji (solid line) 
gives J2 = 1.&6K. 



T = 2K. Good agreement is achieved for J2 = 1.95 K. 
For this value of the exchange couplings, we also plot the- 
oretical magnetization curves at T = 8K and T = 20 K 
in Fig. These curves deviate significantly from the 
data, but one should be aware that T = 8K is already 
beyond the estimated limit of validity T < 0.5 JqS of our 
theoretical approach. 

In Fig. 1111 the uniform susceptibility is plotted in the 
experimental units cm 3 /mol. When all exchange inte- 
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grals have the same order of magnitude we expect a peak 
in the susceptibility for T ~ JqS. Experimentally, the 
peak is at approximately 7 K so that we have Jo ~ 3 K, in 
accordance with the fits of the magnetization curves. For 
a more quantitative comparison we use the following pro- 
cedure. First we substract the temperature-independent 
contribution from the experimental susceptibility in or- 
der to get the correct paramagnetic behavior at high tem- 
peratures. Then we fit the theoretical expression (|3.39|) 
with J 2 = 2 Jx to the full set of data points. Circles in 
Fig. ^2 are experimental data and the solid line is a fit 
with J2 = 1.66 K. The theoretical curve reproduces the 
behavior of the susceptibility very well and it especially 
gives a good estimate of the position and the form of the 
peak. Note that we experimentally observe an increase 
in the susceptibility below T* = 3.0 ± 0.2K. This coin- 
cides with an anomaly in the specific heat. The care- 
ful reader will notice at this point that the estimated 
value of T» is larger than the temperature T — 2 K where 
we obtained the best fit of our calculated magnetization 
curve m(H) to the experimental data shown in Fig. ^| 
Hence, at T = 2 K the system seems to have some kind 
of long range magnetic order, which we have ignored in 
our calculation. However, the precise nature of the order 
and the mechanism responsible for the ordering are not 
known at this point. The fact that a strictly 2D model 
can reasonably well explain the magnetization curve at 
T = 2 K imposes some constraint on possible ordering 
mechanisms. We suspect that dipole-dipole interactions 
play an important role in this temperature range, because 
the long-range nature of the dipole-dipole interaction can 
give rise to spontaneous antiferromagnetic order even in 
2D ml This point deserves further attention, both theo- 
retically and experimentally. 



C. Staggered correlation length in a magnetic field 

The energy gap appearing in Eq. I|3.23|l can be related 
to the staggered correlation length £, as discussed by 
Takahashi. 17 Assuming for simplicity \S V \ = a, we may 
identify 



4ft. 
n 



(4.4) 



In the absence of a uniform field the low temperature be- 
havior of £ has been thoroughly studied by Chakravarty, 
Halperin and Nelson^ Surprisingly, the effect of a uni- 
form field h on £ has so far not been investigated. We 
now analyze the asymptotic behavior of £ at low tempera- 
tures. In two spatial dimensions, the limit T — > also im- 
plies h s — ► 0. Our self-consistency equations l|3.25|l and 
(|3.26l) can then be solved analytically by isolating diver- 
gent contributions to the integrals I(h, h s ) and F(h, h s ) 
originating from gapless modes in the spin-wave spec- 
trum. In the regular part of the integral, the limit T — > 
and h s — > can be taken. For the leading behavior at 
small uniform fields h <C 1 only the singular part of 



I(h, h s ) contributes, and we obtain the self-consistency 
condition 



= n(0) 



I sin£ (K h s ) 
S 



(4.5) 



Here, I slns (h, h s ) is the part of the integral I(h, h s ) that 
diverges for vanishing gaps in the spin- wave dispersions, 
and n(0) = n(h = 0,h s = Q,T = 0). For h < 1, we 
obtain 



S (M.) 



E 



T V u 
JoS 2 

T V. 



J S 871-44 

In I Ah 2 



d 2 k 1 
n 



In 



Ahs 
n 



(4.6) 



to leading logarithmic order. From Eqs. I|4.4J) and H4.5(l 
we then obtain the following result for the self-consistent 
energy gap in a small uniform magnetic field 



A 2 (h) = 



2£(h) 



(2h) 4 {2hf 



(4.7) 



where A = a/2£(0) is the gap for vanishing uniform 
field and the temperature dependence of the zero-field 
staggered correlation length is given by 



oc exp 



2irJ S 2 n{0) lj y 



T 



For a square lattice this yields with Jo 

Uy/V u = 1/4 



oc exp 



2irJS 2 n(0) 
f 



(4.8) 
4J and 

(4.9) 



which is identical to Takahashi's result (see Eq. (27a) 
in Ref. I17B . except that we do not include a spin- wave 
velocity renormalization in our approach. To obtain this 
renormalization, the spin-wave interaction would have to 
be treated on the mean-field level in a fully self-consistent 
way. 

The field dependence of the correlation length for fixed 
temperature is given by Eq. I|4.7() . For h -C A (T), we 
have 



m = c(o) 

whereas for h 3> A (T), we obtain 



1 

1+ 2 



h 



4/, i'W 



(4.10) 



(4.11) 



From Eq. (|4~7|l it is clear that £(/i) > £(0). Thus, the 
correlation length is increased by a small uniform field 
due to reduced quantum fluctuations. 
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The temperature dependence of the correlation length 
for fixed uniform field h can also be extracted from 
Eq. I4.7jl . As long as A (T) > 2h, this temperature 
dependence is still given by Eq. 14. 8|) . When the temper- 
ature is further reduced, Eq. (|4.7|l predicts a crossover 
at Ao(T) « 2h to the following temperature-dependent 
correlation length 

£(h) f A-irJ S 2 n(0) £J y \ 



The additional factor of two in the exponent as com- 
pared to Eq. H4.8fl is due to the fact that at very low 
temperatures the spin-wave mode yields a singular 
contribution, whereas the mode €k+ has a gap 2h which 
is fixed by the external field. In contrast, for h — both 
modes contribute equally, leading to Eq. (|4.8|l . 

The analysis in this section has been carried out for 
h <C 1. For larger fields, there are field dependent pref- 
actors of the first logarithm in Eq. (|4.6|) leading to a field 
dependent renormalization factor Zh in the exponent of 
Eq. i|4.12|l . The field dependence of the correlation length 
at fixed temperature is then no longer determined by the 
singular contributions to the integrals and cannot be ex- 
tracted from the simple analysis presented here. Close to 
the critical field at h = 1 the nature of the divergences 
changes, since the dispersion of the a = — mode becomes 
quadratic. As our mean-field calculation is not suitable 
to describe the true critical behavior in two dimensions, 
we do not discuss this limit in more detail. 

Our approach can also describe a quasi one- 
dimensional anisotropic system, where the exchange cou- 
pling between chains is very weak. The dispersion is then 
almost flat in the transverse direction. The integrals will 
be quasi one-dimensional as long as the maximum varia- 
tion of the dispersion in the transverse direction is smaller 
than the self-consistent gap 4h s /no. In this intermediate 
temperature regime the staggered correlation length be- 
haves as if the system were one-dimensional. At even 
lower temperatures there will be a crossover to the true 
asymptotic two-dimensional behavior. A rough estimate 



for the position of the crossover is a/£ oc £±/a where l\ 
is the eigenvalue of the matrix A defined in Eq. (|3.18|) as- 
sociated with the eigenvector perpendicular to the chain 
direction. 



V. CONCLUSION 

In summary, we have investigated the magnetic 
properties of the new metal-organic quantum magnet 
Mn[Ci H 6 (OH)(COO)] 2 x 2H 2 0. Its layered structure 
contains two-dimensional arrangements of Mn 2+ ions 
that suggest a spin 5 = 5/2 Heisenberg model on a dis- 
torted honeycomb lattice as a minimal model. In order 
to explain measurements of the magnetization M (H, T) 
and the susceptibility x(T), we develop a variant of mod- 
ified spin-wave theory, which can be used to describe fi- 
nite temperature properties of two-dimensional magnets 
in a uniform external magnetic field. A fit of the theo- 
retical results to the experimental curves shows a satis- 
factory agreement for the magnetization at low temper- 
atures where we expect our theoretical approach to be 
valid. The magnetic susceptibility is very well described 
down to temperatures of T > X* w 3K. Both quanti- 
ties are consistently fitted by one parameter J 2 = 2Ji 
to give the exchange coupling J 2 ~ 1.8 K. For tempera- 
tures below T* the uniform susceptibility shows again an 
upturn, which together with an anomaly in the specific 
heat is most likely due to some ordering transition. Pos- 
sible mechanisms for this transition are dipole-dipole in- 
teractions or couplings between the layers, which should 
be included in more refined theoretical models. From 
the experimental point of view nuclear magnetic reso- 
nance or neutron scattering measurements could provide 
a more detailed insight into the nature of the magnetic 
interactions. 

This work was supported by the DFG via Forscher- 
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importance of dipole-dipole interactions. 
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